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COMPUTATIONAL DESIGN OF LOW ASPECT RATIO 
WING-WINGLET CONFIGURATIONS FOR TRANSONIC 
WIND-TUNNEL TESTS 

by 

John M. Kuhlman and Christopher K. Brown 
SUMMARY 

A computational design has been performed for three different low aspect 
ratio wing planforms fitted with nonplanar winglets; one of the three 
planforms has been selected to be constructed as a wind tunnel model for 
testing in the NASA LaRC 7 x 10 high speed wind tunnel. A design point of M 
= 0.8, Cl a 0.3 was selected, for wings of aspect ratio equal to 2.2, and 
leading edge sweep angles of 45* and 50*. Winglet length is 15% of the wing 
semispan, with a cant angle of 15*, and a leading edge sweep of 50*. Winglet 
total area equals 2.25% of the wing reference area. This report summarizes 
the design process and the predicted transonic performance for each 
configuration. 

INTRODUCTION 

Winglets have proven to be effective nonplanar drag reduction devices in 
several applications to high aspect ratio wing planforms typical of transport 
or business jet aircraft. However, recent studies have indicated even larger 
potential benefits may be obtained when winglets are used on low aspect ratio 
configurations such as fighter aircraft (Refs. 1-3). It was found in the 
computational work of Refs. 1-3 that one can obtain the same percentage 
reduction in drag coefficient at the same Cl and ratio of winglet 
length-to-wing span, independent of wing aspect ratio and leading edge 
sweep, even at the transonic design point selected for the current work. 
Since a low aspect ratio wing has a lower lift-to-drag ratio than a high aspect 
ratio wing, then an equal percentage reduction in drag coefficient at equal lift 
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coefficient results in a larger drag force reduction at low-aspect ratio. 

DESIGN PROCEDURE 

The present work has been undertaken to design a low aspect ratio 
wing-winglet wind tunnel model to be constructed and tested in the NASA 
LaRC 7 x 10 high speed tunnel, to confirm the numerical drag reduction 
predictions of Refs. 1-3. Designs have been performed for three different 
wing planforms, all using the same design procedure developed in Refs. 1-2. 
For each wing planform, an optimum wing-alone geometry and a wing-winglet 
geometry have been defined. A linear potential flow theory design code (Refs. 
4, 5) has been used to define wing-winglet and wing-alone camber surfaces 
for minimum induced drag at the selected design point of M = 0.8, Cl a 0.3. 
This design point was chosen as being representative of a cruise condition for 
heavily loaded lightweight fighters at an altitude of 30,000 feet. The design 
code was run at Cl = 0.4, because addition of a fuselage was found to reduce 
the calculated Cl by approximately 0.1. A NACA 64A006 thickness distribution 
has been added to the camber surface, and a cylindrical fuselage having a 
diameter equal to 0.125b, and 5.25b in length has been used. For all 
wing-winglet configurations the wing and winglet geometry have been altered 
in the vicinity of the wing-winglet juncture, to reduce loading and eliminate 
or reduce the strength of any shocks formed in this region. Wing tip airfoil 
camber has been reduced, and geometric incidence has been reduced for the 
outboard 10% of the wing, while winglet toe out has been increased at the 
winglet root, following the procedures which were developed in Refs. 1-3. 
Also, for all current designs an a = 0.8 chordwise loading shape function has 
been utilized in an effort to elliminate any predicted upper surface trailing 
edge boundary layer separation, such as was found for the earlier designs 
which used an a = 1.0 rectangular loading (Ref. 1). This procedure was 
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successful at eliminating predicted upper surface boundary layer separation 
for an aspect ratio 2.20 wing-winglet and wing in Ref. 3. Figs. 1 and 2, taken 
from Ref. 3, summarize these results. Typical pressure coefficient 
distributions for the A = 2.20, A = 45* wing-winglet are shown in Fig. 1 for 
designs using a = 1.0, 0.9, and 0.8 chordwise loading functions at M = 0.8, a = 
0*. Pressure recovery on the upper surface is more gradual as the value of 
a is reduced, but shocks on the winglet are strengthened slightly. As stated 
in Ref. 3, there was no predicted boundary layer separation on the wing for 
the a = 0.8 configuration. As shown in Fig. 2, all three chordwise loading 
functions yielded essentially the same calculated drag polars. The 
wing-winglet geometry for this planform is shown in Fig. 3. 

Performance predictions for the wing-alone and wing-winglet 
configurations versus angle of attack have been obtained using the 
WIBCO-PPW transonic small disturbance code of Refs. 6,7, at M = 0.8. 

Calculated force coefficients, spanloads, and boundary layer separation 
locations on the wing will be presented for all three wing planforms and the 
three corresponding wing-winglet configurations. Also, typical calculated 
pressure coefficient distributions will be presented. 

PLANFORM DESCRIPTIONS 

Two planforms previously studied in Ref. 1 have been used in the present 
design effort. These wing planforms were called cases F and G in Ref. 1. 
Also, a third planform has been studied which is essentially configuration G 
with an unswept trailing edge (called cropped delta G in this study). 
Definition of these three wing planforms is given in Table 1, while Figs. 3-5 
show the resulting wing-winglet design geometries without the fuselage for 
cases F, G, and cropped delta G, respectively. Wing F has a leading edge 
sweep of 45* and A = 2.2, while wing G has A = 50*, A = 2.2. Both of these 
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wings have a taper ratio of 0.2. The cropped delta G has an unswept trailing 
edge, with A = 50*, A = 2.22, and X = .203. 

All 3 wing-winglet configurations have winglet planforms with A s 50* and 
a taper ratio of 0.5. Winglet root chord is 60% of the wing tip chord and 
winglets have been mounted in an aft position. Winglet cant has been fixed at 
15* from the vertical, and all winglets have used a NACA 64A006 thickness 
distribution. These winglet planform choices are similar to those used in refs. 
1-3, and are similar to design recommendations by Whitcomb (Ref. 8) for 
winglets mounted on transport type wings. Winglet total area is 2.25% of the 
wing reference area for configurations F and G, and 2.27% for configuration 
cropped delta G. 

Wing-alone design geometries obtained from the linear design code have 
not been altered. However, in order to obtain successfully converged 
transonic flow predictions for the wing-winglet geometries using the 
WIBCO-PPW code, it was necessary to modify the linear theory camber surfaces 
in the wing-winglet juncture region, as discussed in Refs. 1,2, to reduce 
loading in this region. In addition, for the a = 0.8 chord loading used in the 
present study it was necessary to further reduce loading in the wing-winglet 
juncture as shown in Table 2. 

PRESENTATION OF PERFORMANCE RESULTS 

Predicted performance results to be presented include lift and drag 
coefficients, pitching moment and wing root bending moment coefficients, 
typical pressure coefficients, normalized spanloads, and upper surface wing 
boundary layer separation locations for the wing-alone and the wing-winglet 
configurations for all three wing planforms, generally for -4* < a < 0*. Note 
that all force and moment coefficients presented include only the forces and 
moments on the wing and winglet, but omit those on the fuselage. Viscous 



effects on the winglet are estimated using an empirical skin friction 
correlation. All results for configurations G and cropped delta G have been 
obtained using 150 crude grid iterations, followed by 150 crude-fine grid 
iterations using the interacted Bradshaw strip boundary layer on the wing, at 
a Reynolds number of 3.8 x 10® based upon wing mean aerodynamic chord. 
This is estimated to be a realistic Reynolds number for the wind tunnel model. 
Results for configuration F (partly taken from Ref. 3) have been obtained 
using 100 crude grid iterations, followed by 200 crude-fine grid iterations 
with the interacted strip wing boundary layer at a Reynolds number of 9 x 
10®. Boundary layer transition has been assumed to occur at x/c = 0.05. 

Note that for both the wing-alone and wing-winglet configuration F, no 
converged solutions could be obtained for « > 0.5*, while all G and cropped 
delta G configurations would not converge for a > 0*. The same difficulty was 
encountered in Refs. 1-3 for the previous geometries using an a = 1.0 

rectangular chord loading. However, for the present configurations using a = 
0.8 chord loadings this difficulty in obtaining converged solutions while 
including the viscous boundary layer calculation at higher lift coefficients 
seems to be worsened. Note also that results for a modified cropped delta G 
wing-winglet are presented (Table 2), where further unloading of the 
wing-winglet juncture by increased winglet root toe out and wing tip twist 
was successful at increasing the angle of attack for which converged solutions 
could be obtained up to « = 1* (Cl = 0.3344, versus Cl = 0.2934 at « = 0*). 

The calculated performance prediction results are presented in the 
following figures: (Force coefficients are also presented in Table 3) 

Figure 

Results Numbers 

Configuration F Force Coefficients 6 

Configuration F Spanloads 7 
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Configuration F Boundary Layer Separation 8 

Configuration F Wing-Winglet Cp’s 9-11 

Configuration G Force Coefficients 12 

Configuration G Spanloads 13 

Configuration G Boundary Layer Separation 14 

Configuration G Wing-Winglet C p ’s 15-17 

Configuration G Wing-Alone Cp’s 18-20 

Configuration Cropped Delta G Force Coefficients 21 

Configuration Cropped Delta G Spanloads 22-23 

Configuration Cropped Delta G Boundary Layer Separation 24-25 

Cropped Delta G Wing-Winglet Cp’s 26-28 

Wing of Cropped Delta G Wing-Winglet Cp’s 29-31 

Modified Cropped Delta G Wing-Winglet C p ’s 32-35 

Optimum Wing-Alone Cropped Delta G Cp’s 36-38 


DISCUSSION OF RESULTS 

Force Coefficients 

Predicted lift and moment coefficients for each of the three different basic 
wing planforms all look quite similar, and all vary linearly versus angle of 
attack. Generally, wing-winglet configurations develop slightly less lift at the 
same « than the corresponding optimum wing-alone configuration. This is due 
to the modifications which were required in the wing-winglet juncture to 
reduce loading in this region. Note, however that the effect of adding a 
winglet to a fixed wing geometry may be seen in Fig. 21 by comparing the 
wing-winglet Cl with that of the wing of the wing-winglet design (diamond 
symbols). Addition of a winglet not only reduces Cp somewhat due to the 
thrust on the winglet but also increases Cl by about 5% (by .016 at Cl a .293 
and by .011 at C L 8 .206). 

Drag polars and L/D versus Cl also all look similar, where the drag polars 
appear to be shifted downwards to lower drag levels for the wing-winglet 
configurations relative to the corresponding optimum wing-alone 
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configurations. Predicted percentage reductions in Cp at equal Cl are 
presented in Table 4 for all three wing planforras for Cl between 0.18 and 
0.26. Note that predicted percent reduction in Cp tends to decrease slightly 
as Cl increases, and that these percent reductions are comparable for all 
three wing planforms. 

Pitching moment coefficients about the wing apex are not altered greatly 
for wing-winglet configurations. For example, for the cropped delta G, C m 
is increased 1.5% at Cl = .18 and 1.9% at Cl = .26 for the wing-winglet 
relative to the wing-alone. Wing root bending moment coefficients for 
wing-winglet configurations are increased by about 5-6% relative to the 
corresponding wing-alone case at equal lift. For the cropped delta G, 
increases are 5.4% at Cl = .18 and 6.0% at Cl = .26. These percentage 
increases are consistent with those observed in Refs. 1-3, and are expected to 
be related to the wing structural weight penalty due to the winglet. 

Spanloads 

Predicted spanload distributions have been normalized by Cl, so total area 
under the curve should equal one and be independent of angle of attack. 
Spanload results are shown typically at « = 0*, -2*, -4*, for both the 
wing-alone and the wing-winglet configurations. Spanload shape does not 
change greatly with angle of attack for the wing-alone configurations. 
However, there is a noticeable, consistent trend for all wing-winglet spanloads 
for the normalized loading to be reduced as angle of attack is increased in 
the vicinity of the wing-winglet juncture, with a corresponding increase in 
inboard loading. Loading near the centerline is reduced for all configurations, 
due to the fuselage. This shifts the loading center outboard, and results in 
higher local Mach numbers on the wing upper surface near the wing tip than 
otherwise would be required to develop a given Cl value. Loading is higher 
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than elliptical near the wing tip of the wing-alone configurations. Also, 

loading outboard on the wing is higher than the theoretical optimum for 
wing-winglet configurations, except at the wing tip. Loading is also reduced 
relative to the linear theory theoretical optimum at the winglet root. Similar 
trends were observed in Refs. 2,3. 

Boundary Laver Separation 

Comparison of predicted wing boundary layer separation locations shows 
a great deal of difference between the three wing planforms. None of the 
configurations have any predicted boundary layer separation on the wing 
lower surface for -4* < « < 0*. The configuration F wing-alone and 
wing-winglet results show essentially no predicted upper surface boundary 
layer separation, as first reported in Ref. 3. However, both the G and 
cropped delta G configurations have predicted boundary layer separation on 
the wing upper surface, which tends to worsen as « is reduced. 

This is shown most clearly in Fig. 39 where predicted upper surface 
separation locations at M = 0.8, « = 0* are compared for the three wing-alone 
designs and the four wing-winglet designs. Note that results for the modified 
cropped delta G wing-winglet are for « = -0.5*, because the solution at « = 0* 
experienced difficulties in the boundary layer calculation. Neither the F or G 
wing-alone configurations have any predicted boundary layer separation, while 
the cropped delta G wing has predicted separation near the trailing edge for 
0.23 < v < .34 and .66 < v < .9. No boundary layer separation is observed for 
the configuration F wing-winglet, while both configurations G and cropped 
delta G have predicted boundary layer separation over the entire wing, from 
the wing-body juncture, where separation is predicted at x/c * .985 to the 
vicinity of the wing-winglet juncture where separation is predicted at x/c a 
.93. The modified cropped delta G is somewhat better, but still has predicted 
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upper surface boundary layer separation outboard of v = .76. 

This trend of increasingly worse trailing edge boundary layer separation 
as the wing trailing edge sweep is increased (from -12* for configuration F, to 
-1.2* for G, to 0* for the cropped delta G) is consistent with trends observed 
for shock-induced trailing edge separation on a series of arrow wings at M > 
1 in Refs. 9,10. However, there is no evidence of any trailing edge shock for 
any of the present results. Hence, the predicted boundary layer separation 
locations from the WIBCO-PPW code have been monitored versus the iteration 
count, as summarized in Table 5. Generally, the predicted separation region 
initially grows and then decreases in size as the iteration count increases. 
The predicted separation region may be further reduced in size or even .be 
eliminated with a greater number of iterations. Also, it is expected that the 
reduced number of initial crude grid iterations used for configuration F (100 
versus 150) may have influenced the boundary layer separation prediction, by 
reducing the steepness of any regions of rapid pressure recovery. 

In an effort to obtain an independent measure of the reliability of the 
boundary layer separation prediction for the G and cropped delta G 
configurations, the shock-induced trailing edge separation criterion developed 
by Cunningham, et al,, (Ref. 11) has been used to analyze the airfoil geometry 
of these configurations, as shown in Table 6. Here the incidence angle for 
onset of shock-induced trailing edge separation, «ted» tabulated, calculated 
according to the method developed in Ref. 11. Airfoil geometric parameters 
used in this method include the nondimensional radius of curvature at the 
airfoil upper surface crest, c/R, and the lower surface and camber line slopes 
at the airfoil trailing edge, <5tel a nd tf tec> respectively. Here, if the value of 
«^ e( j is negative, then shock-induced trailing edge separation i3 predicted to 
have already occurred, while «ted > 0 indicates that separation should not 


9 



occur until « > «ted* Comparison of Table 6 and Fig. 39 shows that the 
separation prediction using Ref. 11 shows the same trend as the strip 
boundary layer calculation in WIBCO-PPW. 

Pressure Coefficient Distributions 

Pressure coefficient distributions for all configurations appear quite 
similar to one another at nearly equal Cl values. Also, Cp distributions on 
the wing of each wing-winglet configuration are essentially identical to those 
of the corresponding wing-alone configuration except at the wing tip, where 
the presence of the winglet results in additional loading. Pressure 
distributions at « = 0* all are quite similar to those obtained previously in 
Refs. 1-3. Use of the a = 0.8 chordwise loading function results in more 
gradual pressure recovery on the upper surface near the trailing edge 
relative to results with a = 1.0, as seen previously in Ref. 3. Mid-chord 
shocks are found on the inboard surfaces of the lower half of the winglets for 
all four wing-winglet configurations for « > -2* (Cl > 0.2). Calculated 

pressure distributions change quite significantly versus angle of attack, even 
though there is not much variation in normalized spanload. 

As angle of attack is decreased to « = -4*, pressure suction spikes are 
observed on the lower surfaces of all wing-alone and wing-winglet 
configurations near the leading edge. The level of these suction spikes at « = 
-4* appear to be quite similar for all 3 wing-alone and 4 wing-winglet 
configurations. The development of such leading edge suction spikes is due 
to the relatively small nose radius of the 64A006 thickness distribution utilized 
for the present design geometries. 

Pressure distributions for the two wings analyzed for the cropped delta 
G planform are quite similar. However the wing of the wing-winglet design 
has slightly greater suction on the upper surface near the leading edge. 
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Calculated upper and lower surface pressure distributions on the winglet, 
and on the wing near the tip, are at times observed to cross near the trailing 
edge. This is believed to be due to the frozen streamwise wake modeling 
utilized in the WIBCO-PPW code. 


CONCLUSION 

Predicted transonic flow performance results have been presented for 
seven different low aspect ratio configurations (three wing designs and four 
wing-winglet designs) for a design point of M = 0.8, Cl a 0.3. All 

wing-winglet designs yield the same predicted percent drag reduction relative 
to the corresponding wing-alone design. However, since it is felt that the 
cropped delta G wing planform is most representative of wing planforms for 
current and next generation fighter wings, this will be the configuration 
which will be constructed for the wind tunnel test, even though this planform 
had the worst predicted boundary layer separation characteristics. The 
modified cropped delta G wing-winglet and cropped delta G wing-alone 
geometries will be constructed to fit to a simplified cylindrical fuselage with 
an ogive nose, to allow a fair comparison between the drag of the 
wing-winglet relative to the wing-alone. Predicted drag reductions due to the 
winglet of about 12% at Cl = 0.26, neglecting the fuselage, should correspond 
to about a 6-8% total drag reduction when the fuselage forces are included. 

Since the configurations selected for the wind tunnel test do have some 
predicted boundary layer separation, it is recommended that some redesign of 
both configurations be performed using the automated design method of Smith 
(Ref. 12), which uses the methodology of the airfoil design method of Campbell 
(Ref. 13). In particular, pressure recovery should be made more gradual near 
the trailing edge, to eliminate the predicted trailing edge boundary layer 
separation for both the wing-alone and wing-winglet designs. Also, it may be 
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possible to improve flow in the vicinity of the wing-body juncture. It may 
also be desirable to increase the nose radius slightly to reduce the leading 
edge pressure spikes away from the design point. However, this may be a 
disadvantage when the model is tested at supersonic Mach numbers. 

Note that the uncambered airfoils having significant positive geometric 
incidence which are found at the wing root for all of the present designs 
(Figs. 3-5) are similar to those obtained for higher aspect ratio wings using 
automated optimization methods and transonic analysis codes, as found in Refs. 
14,15. In Reference 14 the starting airfoil geometries included aft-cambered 
supercritical sections at the wing root, but the twisted, uncambered final root 
airfoils reduced the configuration drag. 

Future efforts will be aimed at obtained WIBCO-PPW performance 
predictions for the cropped delta G configurations using the actual fuselage 
geometry, once the fuselage geometry has been finalized. Results will be 
obtained both with the viscous boundary layer calculation, as well as without 
the boundary layer to obtain results at higher Cl values. Also, performance 
predictions will be obtained using the cylindrical TAG grid version of PPW 
developed by Rosen (Ref. 16), and the store carriage code of Ref. 17. This 
code utilizes rotated finite differences to better capture shocks at higher Cl 
and can be run at low supersonic Mach numbers. Finally, construction of a 
low aspect ratio wing-winglet model designed at M = 0.1 will be completed, and 
the configuration will be tested in the WVU low speed wind tunnel at velocities 
of 100-200 ft/sec. 
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Table 1. Wing Planform Definition 


1. 

CASE F 

A = 

2.2 

TR = 

0.2 

SWEEP = 

45* 

2. 

CASE G 

A = 

2.2 

TR = 

0.2 

SWEEP = 

50* 

3. 

CROPPED DELTA G 

A = 

2.22 

TR = 

0.203 

SWEEP = 

50* 


Table 2. Incidence Variation for Wing-Winglets at M = 0.8 
Using a = 0.8 Chord Loading 


Change in Change in 

Wing Tip Incidence Winglet Incidence 



at v = (0.91.0.97,1.0) 

at ( = (0, .42. .80,1.0) 

CASE F 

o 

• 

0 

• 

1 

• 

-5*, -3*, -1*, 0* 

CASE G 

-0.6*, -1.2*, -1.3* 

-3.9*, -2.5*, -0.9*, 0* 

CROPPED DELTA G 

-0.6*, -1.2*, -1.3* 

-3.9*, -2.5*, -0.9*, 0* 

MODIFIED CROPPED 
DELTA G 

-1*, -1.8*, -1.9* 

-4.5*, -3*, -0.9*, 0* 
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Table 3. Calculated Force and Moment Coefficients 


Configuration F 


a 

C L 

C D 

C 

m 

C B 

Configuration 

0.0 

.28616 

.01870 

-.2546 

.14628 

FWING 

-1.0 

.24538 

.01380 

-.2265 

. 12551 

FWING 

-2.0 

.20326 

.00984 

-.1971 

.10390 

FWING 

-3.0 

. 16040 

.00700 

-.1670 

.08178 

FWING 

-4,0 

.11432 

. 00568 

-.1345 

.05765 

FWING 

0 . 0 

.27881 

.01536 

♦ Km *m 

.15118 

FWWLT 

-1.0 

.23507 

.01077 

-.2205 

.12803 

FWWLT 

-2.0 

.19042 

.00740 

-.1882 

.10422 

FWWLT 

-3.0 

.14428 

.00541 

-.1548 

.07949 

FWWLT 

-4.0 

.09667 

.00514 

-.1206 

.05371 

FWWLT 


Configuration G 


a 

C L 

C D 

C 

m 

C B 

Configuration 

0.0 

.29192 

.01978 

-.2908 

.15008 

GWING 

-0.5 

.27466 

.01779 

-.2786 

.14144 

GWING 

-1 . 0 

.25110 

.01504 

-.2580 

.12936 

GWING 

-1 .5 

.22984 

.01303 

-.2408 

.11838 

GWING 

-2.0 

.21174 

.01173 

-.2280 

. 10910 

GWING 

-2.5 

. 18979 

.01014 

-.2099 

.09760 

GWING 

-3.0 

.16652 

.00887 

-.1905 

.08545 

GWING 

-3.5 

. 14347 

.00806 

-.1717 

.07355 

GWING 

-4.0 

. 12070 

.00773 

-.1535 

.06153 

GWING 

0.0 

.29390 

.01770 

-.2937 

.15984 

GWWLT 

-0.5 

.27124 

.01516 

-.2744 

.14766 

GWWLT 

-1.0 


♦ 

♦ 

♦ 

GWWLT 

-1 . 5 

.22679 

.01102 

-.2373 

. 12376 

GWWLT 

-2.0 

.20481 

► 00954 

-.2194 

.11196 

GWWLT 

-2.5 

. 18225 

.00831 

-.2008 

.09979 

GWWLT 

-3.0 

. 15923 

.00745 

-.1818 

.08705 

GWWLT 

-3.5 

. 13388 

.00685 

-.1601 

.07361 

GWWLT 

-4.0 

.11043 

.00693 

-.1412 

.06060 

GWWLT 
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, Table 3. (Concluded) - Calculated Force and Moment Coefficients 


Configuration Cropped Delta G 


a 

C L 

C D 

C 

m 

C B 

Configuration 

0,0 

.2927? 

.01971 

-.2907 

. 15044 

CUNG0PT 

-0.5 

.27364 

.01746 

-.2761 

. 14097 

CUNGOPT 

-1.0 

.25238 

.01511 

-.2587 

. 1299? 

GWNG0PT 

-1.5 

.23003 

.01295 

-.2402 

.11848 

CUNGOPT 

-2.0 

.21162 

.01165 

-.2270 

.10920 

CUNGOPT 

-2.5 

. 18708 

.00981 

-.2058 

.09622 

CUNGOPT 

-3.0 

. 16658 

.00890 

-.1902 

.08556 

CUNGOPT 

-3.5 

.1429? 

.00808 

-.1708 

.07332 

CUNGOPT 

-4.0 

.11959 

.00777 

-.1519 

.06092 

CUNGOPT 

0.0 

.2776? 

.01844 

-.2733 

. 14335 

CUNGNOT 

-0.5 

.25815 

.01616 

-.2582 

. 13363 

CUNGNOT 

-1.0 

.23633 

.01379 

-.2399 

. 12250 

CUNGNOT 

-1.5 

.21561 

.01194 

-.2235 

.11216 

CUNGNOT 

-2.0 

. 19506 

.01035 

-.2074 

. 10147 

CUNGNOT 

-2.5 

. 17402 

.00900 

-.1908 

.09088 

CUNGNOT 

-3 . 0 

. 15088 

.00778 

-.1715 

.07826 

CUNGNOT 

-3.5 

.12788 

.00707 

-.1528 

.06670 

CUNGNOT 

-4.0 

♦ 10475 

.00676 

-.1342 

.05445 

CUNGNOT 

0.0 

.29344 

.01730 

-.2949 

♦ 1595? 

CULTOLD 

-0.5 

.27181 

.01488 

-.2769 

.14814 

CULTOLD 

-1 . 0 

♦ 

% 


♦ 

CULTOLD 

-1.5 

.22862 

.01107 

-.2415 

.12477 

CULTOLD 

-2.0 

.20594 

.00944 

-.2224 

.11267 

CULTOLD 

-2.5 

. 18278 

.00821 

-.2031 

. 10011 

CULTOLD 

-3.0 

. 16091 

.00754 

-.1857 

,08840 

CULTOLD 

-3.5 

,13517 

.00688 

-.1632 

.07429 

CULTOLD 

-4.0 

.11080 

.00692 

-.1430 

.06096 

CULTOLD 

1.0 

.33443 

.02268 

-.3286 

. 18121 

CULTNEU 

-0.5 

. 26995 

.01472 

-.2748 

. 14687 

CULTNEU 

-1.0 

.24789 

.01259 

-.2564 

. 13492 

CULTNEU 

-1.5 

.22561 

.01084 

-.2379 

. 12261 

CULTNEU 

-2.0 

,20356 

.00937 

-.2197 

.11107 

CULTNEU 

-2.5 

. 18022 

, 00831 

-.2006 

. 09806 

CULTNEU 

-3.0 

. 15831 

.00754 

-.1826 

.0865? 

CULTNEU 

-3.5 

. 13301 

.00701 

-.1608 

.07264 

CULTNEU 

-4.0 

.1083? 

.00712 

-.1403 

.05918 

CULTNEU 
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Table 4. WIBCO-PPW Predicted Percentage Drag Reductions Due to Winglets 
at M = 0.8 Using a = 0.8 Chord Loading 



C L = 0.18 

C L = 0.22 

C L = 0.26 


CASE F 

18% 

16% 

12.7% 

CASE G 

14.8% 

14.6% 

12.7% 

CROPPED DELTA G 

14.7% 

15.4% 

13.3% 
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Table 5. Predicted Upper Surface Boundary Layer Separation Locations 
Versus Iterations for Modified Cropped Delta G Wing-Winglet 
at M = 0.8, « = -0.5* (150 crude grid iterations) 

x/c for Boundary Layer Separation 


V 

154 its 

214 its 

254 its 

314 its 

354 its 

414 its 

.145 

.985 






.195 

.983 






.245 

.983 






.295 

.979 






.347 

.980 






.400 

.978 






.455 

.972 






.511 

.969 






.570 

.967 


.999 



.996 

.631 

.965 






.695 

.963 


.972 




.763 

.960 


.957 

.997 

.978 

.990 

.836 

.955 

.979 

.955 

.976 

.948 

.979 

.914 

.956 

.952 

.951 

.951 

.947 

.957 
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Table 5. (Concluded) - Predicted Upper Surface Boundary Layer Separation 
Locations Versus Iterations for Cropped Delta G Optimum Wing- 
Alone at M = 0.8, a = -0.5* (150 crude grid iterations) 

x/c for Boundary Layer Separation 


V 

154 its 

214 its 

254 its 

314 its 

354 its 

414 its 

450 its 

.139 

.984 







.186 

.982 







.234 

.980 







.282 

.975 







.332 

.976 






.996 

.382 

.973 



.995 


.994 


.435 

.968 



.993 


.990 


.488 

.967 



.990 


.996 


. 544 

.962 



.985 

.988 



.603 

.955 


.965 

.990 

.996 

.990 


.665 

.951 


.959 

.990 

.995 

.987 


.730 

.945 

.997 


.990 

.994 

.985 

1.000 

.799 

.942 

.992 


.987 


.991 

.993 

.874 

.940 

.988 

.996 


.997 

.992 

.994 

.956 

.951 




.998 
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Table 6. Shock-Induced Trailing Edge Boundary Layer Separation 
Prediction Using Method of Ref. 11 


Configuration 

( x/ c ) crest 

(c/R) 

<5tel(rad) 

<5tec( rad ) 

«t< rad ) 

«ted 

Crop Delta G 
w-wlt 

0.4 

.4584 

.02761 

.09299 

.0494 

-3.509* 

Crop Delta G 
wing 

0.4 

.4812 

.04157 

.10671 

.0512 

-4.309* 

G w-wlt 

0.4 

.4583 

.02770 

.09307 

.0502 

-3.555* 

G wing 

0.4 

.4817 

.04183 

.10695 

.0522 

-4.380* 

F w-wlt 

0.4 

.4404 

.02522 

.09058 

.0534 

-1.758* 

F wing 

0.4 

.4676 

.03887 

.10394 

.0578 

-2.674* 


21 




a = 1.0, M = 0.8, « = 0* 

a = 0.9, M = 0.8, « = 0* 

a = 0.8, M = 0.8, « = 0* 



<• = 0.650 



Fig. 1 


Calculated wing-winglet C 
0.8, Cl * 0.27, for various 


p distributions for configuration F at M 
chord loading shape functions, from ref. 


3. 
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-0 

Wing-Alone, a = 1.0 

6- 

A 

Wing-Winglet, 

a = 1.0 

A 

Wing-Alone, a = 0.9 


Wing-Winglet, 

a = 0.9 

o 

Wing-Alone, a = 0.8 


4 

Wing-Winglet, 

a = 0.8 




Fig. 2 Predicted performance of wing-alone and wing-winglet configurations 
F at M = 0.8, for various chord loading shape functions, from ref. 3. 
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Fig. 6. (Concluded) - Predicted performance of wing-alone and wing-winglet 
configurations F at M = 0.8; drag polar and L/D - Cl« 
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-3.2 

-2.4 

- 1.6 

C -0.8 
P 

0 

0.8 

1.6 



Fig. 9. Calculated wing-winglet C p distributions for configuration F at M 
0 . 8 , « = 0 *. 
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Fig. 10. Calculated wing-winglet C p distributions for configuration F at M = 
0.8, a = -1.5*. 
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ETA WINS 


DIAMOND: A-1 
CIRCLE: A-0 


TRIANGLE: A— 2 
SQUARE: * 



Fig. 13. Calculated normalized spanloads for wing-alone and wing-winglet 
configurations G at M = 0.8. 
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Fig. 15. Calculated wing-winglet Cp distributions for configuration G at M 
0 . 8 , « = 0 *. 


37 




Fig. 16. Calculated wing-winglet Cp distributions for configuration G at M 
0.8, « = -1.5*. 
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L Wing Upper Surface 





Fig. 18. Calculated wing-alone Cp distributions for configuration G at M = 0.8, 

« = o*. 
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Fig. 19. Calculated wing-alone C p distributions for configuration G at M = 0.8, 
« = - 2 *. 
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Fig. 20. Calculated wing-alone C p distributions for configuration G at M = 0.8, 

a — *4 . 
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ANGLE Of ATTACK (DECREES) 




Fig. 21. Predicted performance of wing-alone and wing-winglet configurations 
cropped delta G at M = 0.8; Cl - a > C m - Cl» Cjj - Cl* 
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LIFT COEFFICIENT 

Fig. 21. (Concluded) - Predicted performance of wing-alone and wing-winglet 
configurations cropped delta G at M = 0.8; drag polar and L/D - Cl. 
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ETA 


Fig. 22. Calculated normalized spanloads for wing-alone and wing-winglet 
configurations cropped delta G, wing of wing-winglet and 
wing-winglet at M = 0.8. 
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DIAMOND: A-1 TRIANGLE: A— 2 

CIRCLE: A-0 SQUARE: A— 4 



Fig. 23. Calculated normalized spanloads for wing-alone and wing-winglet 
configurations cropped delta G, optimum wing and modified 
wing-winglet at M = 0.8. 
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Fig. 24. Predicted upper surface boundary layer separation locations for 
wing-alone and wing-winglet configurations cropped delta G, wing of 
wing-winglet and wing-winglet at M = 0.8. 


47 




0 - 0 - 0 0 0-0-0 - 2.3 A-A A -3 


O O O -3.5 B-Q-Q -4 



■INC SEMI-SPAN LOCATION 


Fig. 25. Predicted upper surface boundary layer separation locations for 
wing-alone and wing-winglet configurations cropped delta G, optimum 
wing and modified wing-winglet at M = 0.8. 
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Fig. 26. Calculated wing-winglet Cp distributions for configuration cropped 
delta G at M = 0.8, « = 0*. 
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Fig. 27. Calculated wing-winglet Cp distributions for configuration cropped 
delta G at M = 0.8, « = -2*. 
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Fig. 28. Calculated wing-winglet C p distributions for configuration cropped 
delta G at M = 0.8, « = -4*. 
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Fig. 29. Calculated wing-alone C p distributions for wing of cropped delta G 
wing-winglet at M = 0.8, « = 0*. 
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Fig. 30. Calculated wing-alone C p distributions for wing of cropped delta G 
wing-winglet at M = 0.8, a = -2*. 
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Fig. 31. Calculated wing-alone Cp distributions for wing of cropped delta G 
wing-winglet at M = 0.8, a = -4*. 
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Fig. 32. Calculated wing-winglet C p 
cropped delta G at M = 0.8, 



distributions for configuration modified 
= + 1 \ 
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Fig. 33. Calculated wing-winglet C p distributions for configuration modified 
cropped delta G at M = 0.8, « = -0.5*. 
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Fig. 34. 


Calculated wing-winglet C p distributions for configuration modified 
cropped delta G at M = 0.8, « = -2*. 




Fig. 35. Calculated wing-winglet C p distributions for configuration modified 
cropped delta G at M = 0.8, ot = -4*. 







Fig. 36. Calculated wing-alone Cp distributions for configuration cropped delta 
G at M = 0.8, a = 0*. 
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Fig. 37. Calculated wing-alone C p distributions for configuration cropped delta 
G at M = 0.8, « = -2*. 
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Fig. 38. Calculated wing-alone Cp distributions for configuration cropped delta 
G at M = 0.8, « = -4*. 
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